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Pion interferometry ("HBT") measurements relative to the reaction 
plane provide an estimate of the transverse source anisotropy at freeze-out, 
which probes the system dynamics and evolution duration. Measurements 
by the STAR Collaboration indicate that the source is extended increas- 
ingly out-of-plane with increasing impact parameter, suggesting a short 
evolution duration roughly consistent with estimates based on azimuthally- 
integrated HBT measurements. 
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1. Introduction 

Two-particle intensity interferometry (Hanbury-Brown-Twiss, or HBT) 
analyses of relativistic heavy ion collisions represent the most direct exper- 
imental probes of the space-time structure of the source generated in such 
collisions 1 . At the Relativistic Heavy Ion Collider (RHIC), identical-pion 
HBT studies in Au+Au collisions [H El |1] yielded an apparent source size 
quantitatively consistent with measurements at lower energies, in contrast 
to predictions of larger sources based on QGP formation 5 . In addition, 
hydrodynamical models, successful at RHIC in describing momentum-space 
quantities such as transverse momentum spectra and elliptic flow have 
failed to reproduce the small HBT radii 0. This "HBT puzzle" 0EI may 
arise because the system's lifetime is shorter than predicted by models. 

In non-central collisions, the initial anisotropic collision geometry gener- 
ates greater transverse pressure gradients in the reaction plane than per- 
pendicular to it. This leads to preferential in-plane expansion (elliptic 
flow) 1101 I1H I12j which diminishes the initial spatial anisotropy. The 
final (freeze-out) source shape should be sensitive to the evolution of the 
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pressure gradients and the system lifetime; a long-lived source would be 
less out-of-plane extended and perhaps in-plane extended. Hydrodynamic 
calculations .13] predict a strong sensitivity of the HBT parameters to the 
early source conditions and show that, while the source may still be out- 
of-plane extended after hydrodynamic evolution, a subsequent rescattering 
phase |14j tends to make the source in-plane. Thus, the freeze-out source 
shape might discriminate between scenarios of the system's evolution. 

Here, we present results of a measurement of azimuthally-sensitive HBT 
in Au+Au collisions at ^/snn = 200 GeV. We discuss how these results may 
be used to extract, in a model-dependent way, the geometric anisotropy of 
the source at freeze-out. Finally, we use this extracted geometry to estimate 
the evolution timescale of the system. 

2. Measurement 

The measurements were made using the STAR detector j!5| at RHIC. 
Particle trajectories and momenta were reconstructed using a Time Pro- 
jection Chamber (TPC) with full azimuthal coverage, located inside a 0.5 
Tesla solenoidal magnet. Au+Au event samples of 0-5%, 5-10%, 10-20%, 
20-30% and 30-80% of the total cross hadronic section (based on charged 
particle multiplicity) are presented here. The second-order event plane an- 
gle ^2 |16| for each event was determined from the weighted sum of primary 
charged-particle transverse momenta (fir) |17|. Within the resolution which 
we determine from the random subevent method ^HJj ^2 ~ ^rp (true re- 
action plane angle) or ^ ~ ^ T p + n", i-e. the direction of the impact 
parameter vector is determined up to a sign [TS|, and the measured event 
plane is roughly coplanar with the true reaction plane 18 . 

Pions were selected according to their specific energy loss (dE/dx) in 
the TPC in the rapidity range |y| < 0.5. Pair-wise cuts to remove two- 
track merging and single-track splitting effects were applied to signal and 
background distributions [2J. 

Pairs of like-sign pions were placed into bins of = ^ pa i r — ^2> where 
^pair is the azimuthal angle of the total pair momentum (k = pi + P2)- 
Because we use the 2 nd -order event plane, is only defined in the range 
(0, 7r). For each bin, a three-dimensional correlation function is constructed 
in the "out-side-long" decomposition ^5] of the relative momentum q. 

Finite reaction plane resolution and finite width of the <I> bins has the ef- 
fect of reducing the measured oscillation amplitudes. A model-independent 
correction procedure |2U| is applied to each q-bin in the numerator and 
denominator of each correlation function, the overall effect of which is to in- 
crease the amplitude of the oscillations of the HBT radii vs. $ (~ 10 — 30%). 

To account for the effect of final-state Coulomb interactions, we adopt 
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the approach first proposed by Bowler an d Sinyukov |22| and recently 
advocated by the CERES collaboration j^Sj- We fit each experimental cor- 
relation function to the form: 

C(q, $) = N ■ [(1 - A) • 1 + A • K(q)(l + G(q, <&))], (1) 

where the (1 — A) and A terms account for the non-participating and partic- 
ipating fractions of pairs, respectively, K{Q im ) is the square of the Coulomb 
wave- function, TV" is the normalization, and G(q, <!>) is the Gaussian form |19j : 

G(q, $>) = e -£ R l^-£ B2 s(®)-fi R2 i^)-Ws R ls^\ (2) 

-R? are the squared HBT radii, where the l,s,o subscripts indicate the long 
(parallel to beam), side (perpendicular to beam and total pair momentum) 
and out (perpendicular to qi and q s ) decomposition of q. 



3. Results, Fourier Coefficients and Estimating Eccentricity 

Figure ^ shows the squared HBT radii as a function of $, for three 
centrality classes and 0.15 < kx < 0.6 GeV/c. Curves correspond to a 
Fourier decomposition of the squared-radii oscillations according to |201 124j 

R 2 ( u x _ / <flJ(kr,$)cos(n$)) (/i = o,a,Z) 

TJ - | (#2 (fc T) $ ) sin(n$)) (/x = os) W 

As expected j2], the th -order Fourier coefficients (FCs), i.e the mean value 
of the radii, indicate larger source sizes for more central collisions. We 
verified that the th -order FC for R^, i? 2 and Rf correspond to HBT radii 
obtained in an azimuthally-integrated analysis. Higher-order (n > 2) FCs 
were found to be negligible. 

HBT radii directly measure homogeneity regions ^HH], which constitute 
only part of the "whole" source; hence, connections to the overall source 
geometry inevitably rely on a model. The so-called blast-wave model |27j 
has been relatively successful in describing soft physics observables (spectra, 
elliptic flow, azimuthally-integrated HBT) at RHIC. In the context of a 
generalized blast- wave, with spatial and flow anisotropies, a recent study [27] 
finds that the transverse source shape is related to Rg(^) as 

= (g " g) 1 ^2 1 1 <2 ( ) 



Where R x and R y are the source length scales in and out of the reaction 
plane, respectively. The approximation becomes exact in the case of van- 
ishing transverse flow |28l I25j . In Figure [21 are the th and 2 nd -order FCs, 
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Fig. 1. (Above) Squared HBT radii mea- 
sured relative to the reaction plane, for 
three selections of event centrality. 

Fig. 2. (At right) Blast-wave [27] calcu- 
lations of Fourier coefficients for the os- 
cillations of squared HBT radii, for three 
sources of varying transverse eccentricity. 
Arrows indicate approximation of Eq. ^ 
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as defined in Equation for three sources as calculated in the blast-wave. 
The sources have the same temperature, radial flow, timescales, and average 
size (all adjusted roughly to the values which best describe RHIC data), but 
differing e. Changing from a round (e = 0) to a deformed source increases 
the oscillations of the transverse radii (and, very slightly- note the scale) 
the longitudinal radius. Equation @] arises from the observation that 2e, 
indicated by small arrows on the right panels, approximately reproduces 
the relative oscillations in the transverse radii. Therefore, relying on the 
blast-wave model for context, we may estimate the freeze-out source shape. 

Figure shows the measured FCs, as a function of centrality, for three 
bins in kx, quantifying the increasing size (left panels) and decreasing 
anisotropy (right panels) as the collision becomes more central. In Figure 0J 
we compare the initial anisotropy, e; n itiai, calculated with a Glauber model 
outlined in reference [25], with efi na i, estimated by applying Equation to 
the data of Figure [3] We note that efi na i has a meaningful sign, and its 
positive values indicate that the timescales and/or flow velocities of these 
collisions are not sufficiently large to fully quench the initial out-of-plane 
geometric anisotropy of the system. 



LisaISMD03 printed on February 8, 2008 



■5 



40 : 
,30; 
20: 
10: 
30 1 



CM <0 r 



10- 
50 - 

40: 

^-30: 

20: 

10: 



• ★ 

★ 



★ ★ 



. 0.15 <k t < 0.25 
* 0.25 <k,< 0.35 
0.35 < k, < 0.60 



tilt 



*.05' 

•o.i 

0.1 

u 

H0.05 1 " 



30 



30 



0.05 



-0.05 



0.1 o^> 



:0.05 33 

0) w 



100 200 300 400 100 200 300 400 

# participants 



Fig. 3. Fourier coefficients for the oscillations of squared HBT radii, for 3 fcr bins, 
versus centrality (number of participating nuclcons) measured by STAR. Left pan- 
els: means (0 th -order FC) of oscillations; right panels: relative 2 nd -order oscillation 
amplitudes. See text for details. 



4. A connection to timescales 

The data in Figure should provide useful constraints on timescales and 
evolution in transport models of heavy ion collisions. As a crude example, 
here we take the model-dependent data in Figure 0] one step further. Mo- 
tivated by hydrodynamical calculations [3U], we assume that the transverse 
flow strength grows roughly linearly with time, and anisotropies in the flow 
field set in very early. Since there is no collective flow initially, the flow ve- 
locity and source extensions in (x) and out of (y) the reaction plane evolve 
with time as 

Here, /3 FO is the flow velocity (at the edge of the source) at freeze-out, and 
To is the evolution time of the system- from initial overlap until freeze-out. 
From Equation |SJ we may calculate e(t) in this toy model. 
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Fig. 4. Final (from HBT) ver- 
sus initial (from Glauber) transverse 
source eccentricity, estimated from 
the data in Fig. and Eq. H The 
line represents e in itiai = efinai- 
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Fig. 5. Toy-model calculations for 
the evolution of the transverse ec- 
centricity as a function of time. See 
text for details. 



For collisions between 10-20% of the cross-section (3 rd centrality bin in 
Figures El and 0} , the estimated initial and final anisotropics are indicated 
on Figure El Estimates for were taken from blast-wave fits |27j to mid- 
central collisions at ^/snn = 130 GeV. In order to evolve e, we must assume 
a value for tq. Blast-wave fits |27j to azimuthally-integrated longitudinal 
HBT radii suggest a rather short evolution time of ~ 9 fm/c. Assuming 
this value, we arrive at a consistent picture: as shown by the solid line in 
Figure [51 at t = tq = 9 fm/c, e has reached its measured value. On the 
other hand, using a significantly different value for To yields inconsistent 
results, as shown by the dashed and dotted lines. Thus, our measurements, 
in the context of this admittedly crude model, seem to provide corroborating 
evidence in the direction of evolution timescales on the order of 9 fm/c. 



5. Summary 

We have presented an analysis of azimuthally-sensitive pion interferom- 
etry for Au+Au collisions at RHIC. We observe clear 2 nd -order oscillations 
of the transverse HBT radii, and extract Fourier coefficients. Geometrically, 
as the impact parameter decreases, the th -order FCs increase, indicating 
larger sources, and the relative oscillations decrease, suggesting rounder 
ones. More quantitative statements rely on model interpretations. 

In the blast-wave model, the relative oscillations provide a good estimate 
of the source ellipticity at freeze-out. In this context, we find that the source 
at all centralities is extended out of the reaction plane (though not as much 
as the initial source), suggesting rather short evolution durations. 

This data should provide strong constraints on dynamic models of heavy 



LisaISMD03 printed on February 8, 2008 



7 



ion collisions, particularly as regards timescales and pressure anisotropies. 
To illustrate the point, we performed a toy model calculation which turned 
out to favor the rather short evolution timescales extracted via other means 
previously. 
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